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I. Introduction 



Various types of correspondence between the Kepler- Coulomb and the isotropic- 
oscillator systems have been extensively investigated since the influential work of Levi- 
Civita early this century IQ. Among the correspondences of interest are mappings 
that can be constructed between the radial equations of the quantum systems. This 
subject was initiated over 50 years ago in a paper by Schrodinger addressing the 
solution of eigenvalue problems by factorization. Schrodinger discovered a connection 
between the radial equation of the three-dimensional quantum Coulomb problem 
and the radial equation of a /^-dimensional quantum harmonic oscillator. Using a 
quadratic transformation in the radial coordinate, he showed that the mapping images 
all the states in the three-dimensional discrete Coulomb spectrum only for oscillators 
with D = 2 OT 4. 

Schrodinger's idea was subsequently rediscovered or investigated by a number of 
authors 0, |, ^, |, 0. An extension relating the radial equations of the d-dimensional 
Coulomb system and the D-dimensional oscillator for the special case of even dimen- 
sions d = D was given in ref . p| . A more general mapping for arbitrary d and even D 
that involves a free parameter was presented in ref. , along with the corresponding 
mappings to the supersymmetric partners of these systems. All these correspondences 
involve oscillators in even dimensions, and they incorporate constraints on the allowed 
range of angular momenta. It is possible in general to map all the states of the d- 
dimensional Coulomb system into half the states of a D-dimensional oscillator, where 
d is greater than one and D must be even. 

Recently, it has been proposed that some restrictions on the dimensions or angular 
momenta can be removed with the introduction of suitable analytical deformations 
called quantum defects in one or both systems ||10|. The motivation for this de- 



rives from the use of supersymmetric quantum mechanics |Tl|, |12[ in the context of 
atomic physics [0], where supersymmetry-based quantum-defect theory (SQDT) Jl^ 
provides an explicit example with direct physical relevance. 

One goal of the present paper is to investigate the issue of relaxing the dimensional 
constraints on the radial correspondences via the introduction of analytical quantum 
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defects. The treatment incorporates not only the Coulomb and oscillator systems 
but also their supersymmetric partners. We show that with a suitable choice of 
defect it is indeed possible to remove restrictions on the mappings. For instance, 
among the examples discussed below is a generalized mapping taking any state in the 
three-dimensional radial Coulomb problem into a state in an analytically modified 
three-dimensional radial oscillator. We also briefly consider the continuum states of 
the two radial systems. 

The focus of this work is the set of radial correspondences as summarized above. 
We do not address here the different issue of obtaining surjective mappings between 
the full D-dimensional oscillator and the full d-dimensional Coulomb systems. This 
interesting question has been addressed by a number of authors, originating with the 
parabolic-coordinate transformation of Levi-Civita that relates D = 2 to d = 2 



and with the mapping of Kustaanheimo and Stiefel |]T5| in their work on celestial 
mechanics that relates D = 4 to d = 3. The latter transformation in particular has 



been much investigated in the quantum context |]T6|, |T8|, |T9[, and in recent years 
extensions connecting D = 8 and d = 5 have been studied ^]. While more 
complete than the purely radial mappings, all these surjective correspondences are 
restricted to a narrow range of dimensions. 

The organization of this paper is as follows. Section II consists primarily of back- 
ground on the supersymmetric radial Coulomb and oscillator systems in arbitrary 
dimensions and the known correspondences between them. It contains key equations 
needed in the subsequent sections and provides a perspective useful for our purposes. 
In section III, we introduce SQDT for the radial Coulomb and oscillator systems in 
arbitrary dimensions, and we define mappings relating these systems to the super- 
symmetric sectors of the associated zero-defect cases. The general correspondence 
between different SQDT for the Coulomb and oscillator systems in arbitrary dimen- 
sions is established in section IV. This permits, for instance, the entire set of Coulomb 
radial states to be injected into a subset of the oscillator radial states for any dimen- 
sions, including odd oscillator dimensions. Section V provides a short discussion of 
some results arising for the continuum Coulomb states. We summarize in section VI. 

To distinguish comparable quantities in the two systems, we adopt the convention 
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that lower-case letters are used for Coulomb-system variables while upper-case letters 
are used for oscillator-system variables. An exception is made in denoting energies, 
for which the symbol E with various sub- and superscripts is used in the Coulomb 
system while F is used in the oscillator system. 

II. Preliminaries 

This section establishes our conventions and presents some preliminary material 
and results. Section IIA begins with definitions and solutions for the radial Coulomb 
problem in arbitrary dimensions, while section IIB similarly treats the harmonic- 
oscillator radial problem. A one-parameter mapping between these systems is pre- 
sented in section IIC. Key equations for supersymmetric quantum mechanics are 
given in section IID. The supersymmetric counterpart of the results in section IIA is 
discussed in section HE, while that of sections IIB and IIC is covered in section IIP. 

IIA. Coulomb Bound States in d Dimensions 

The quantum Kepler- Coulomb system in d dimensions is governed by the hamil- 
tonian 

where is the reduced mass, k is the force constant, and r is the usual radial variable. 
To avoid normalization issues, we assume d > 1. The associated radial equation is 
obtained from the Schrodinger equation hip = Eip by separating the wave function 
ip in generalized polar coordinates, ip{r, 6*1, ... , 9d-i) = s{r-)6{6i, . . . , 9d-i)- Normal- 
izable solutions to the ensuing radial equation are found for discrete eigenenergies 
given by 

where Eq = 2^k'^ /%^, 7 = ^(c? — 3), and n is the principal quantum number taking 
values n = /-|-l,/-|-2, with / the angular-momentum quantum number arising 
from the separation of variables. For c? > 3, / = 0, 1, 2, . . . as usual. When d = 2, the 
angular-momentum quantum number takes the values 0, ±1, ±2, In this case, we 
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define the symbol I to represent the modulus of the angular momentum. Note that 
Eq is the magnitude of the ground-state energy in the lowest dimension considered, 
d = 2. 

To simplify the equations that follow, we introduce a dimensionless radial variable 
y = r/ro, where tq = /i^/2/t/i. It is also convenient for later considerations involving 
supersymmetry to work with a scaled radial function w{y) = y^~^^s{r = Toy), which 
effectively removes the first-order derivative appearing in the radial equation for s. 
The radial equation becomes 

+ -ET } Wd,n,l{y) = ^ ■ (3) 



I dy^ y2 y Eq 

The eigensolutions involve Sonine-Laguerre polynomials and are given by 

wa,nAy) = c,ni y'^'^' exp(^) Lftr ^) (--^) , (4) 
with normalization 
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2rg(n + 7)2'+'i+ir(n + l + d- 2) 
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(5) 



IIB. Oscillator Bound States in D Dimensions 

The quantum hamiltonian for the isotropic harmonic oscillator in D dimensions, 
L> > 1, is 

H--^V' + \Mn^R' , (6) 

where M is the oscillator mass, Q, is the frequency, and R is the usual radial vari- 
able. Separating variables in generalized polar coordinates as before produces a radial 
equation that has normalizable solutions S{R) for energy eigenvalues given by 

F^,r = i^o(2iV + 2r + 3) , (7) 

where T — ^{D — ?,), Fq — ^UQ, is the ground-state energy for the lowest dimension 
D = \, and N is the principal quantum number taking values N = L, L-\-2, L-\-A, . . 
with L the quantized angular momentum arising from the separation of variables. For 
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D > 3, L = 0, 1, 2, . . . as usual. For D = 2, the angular momentum ranges over 0, ±1, 
±2, . . ., and we define L to be its modulus. For D = 1, the only possibilities are L = 
and 1. The corresponding angular variable has two discrete values, distinguishing the 
two orientations of the radial vector. These two cases represent distinct radial systems 
for the one-dimensional oscillator, as opposed to the full one- dimensional oscillator 
with configuration space including both positive and negative coordinate values. For 
convenience in what follows, we define the parity of the radial wave functions as even 
if L = and odd if L = 1. 

1 

Defining the dimensionless variable Y = R/Rq with Rq = {h/MQ)'^ and intro- 
ducing for later convenience the scaled radial function = Y^~^^S{R = RqY), 
the radial equation becomes 

^ (L + r)(L + r + i) _^ _ F 1 j^i^^i _ ^ 



( dY^ F2 FoJ 

The eigenfunctions are 

Wn,NAy) = CDNLY'^^^^'exp{-lY')Lf''V'^\Y') , (9) 

2 2 

with normalization 
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(10) 



With our definitions for L above, these expressions hold for all integral D > 1. Note 
that the D = 1 normalization coefficients differ from the canonical ones by a factor 
of V2 because the above construction produces a normalization on the half line only. 

lie. Mappings between the Coulomb and Oscillator Problems 

The wave function (^ can be mapped to the wave function (H) through the 
quadratic transformation 

(n + 7) 

This correspondence also interconnects the differential equations (^) and (^. The 
explicit relation between eigensolutions is 

WD,NAy) = KdnxY-'^ Wd,n,i ({n + ^)Y^) , (12) 
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where the quantity K^nx maintains the normahzation of the wave functions and is 
given by 

(2n + - 3) 

The quantity A provides an extra degree of freedom in the mapping. 

The ensuing relationships among the dimensionahties and the quantum numbers 
of the two systems are 

D = 2d-2-2X , 
N = 2n-2 + X , 

L = 21 + \ . (14) 

The last of these equations constrains A to be integral. It then follows from the first 
equation that this mapping has image only in the oscillators of even dimension D. 

For given angular momenta / and L, the relation (jM]) between the principal quan- 
tum numbers ensures that the stack n > / + 1 of Coulomb states is in one-to-one 
correspondence with the stack N > L of oscillator states, with ground states coincid- 
ing. This relation between and n determines a condition relating the energies E 
and F of the two systems: 



NT 



. (15) 



The factor of two can be viewed as originating from the scaling of relative to n in 
the second equation in (0). Absorbing it in the definitions of Eq or Fq would change 
equations @ or (||). 

Condition (|1^ shows that successive Coulomb angular momenta / map to every 
second oscillator angular momentum L. The entire set of radial states \n, I) of the 
(i-dimensional Coulomb system can be mapped into a subset of the states |A^, L) of 
the D-dimensional oscillator provided D satisfies 2 < D <2d — 2. For even or odd A, 
the mapping is then an isomorphism to even or odd L, respectively. For given d, the 
allowed values of the pair (D, A) characterizing this mapping are distinct: (2, d — 2), 
{4,d — 3), . . ., {2d — 4, 1), (2d — 2, 0). We recover in this way Schrodinger's result 
that all states of the three-dimensional Coulomb system can be mapped only into 
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oscillators of dimension two or four. Note that if instead / is taken to be fixed, so that 
only a subset of states is imaged, then the allowed range of D is^2 < D < 2d— 2+41. 

IID. Supersymmetric Quantum Mechanics 

For the purposes of the present paper, only a few of the basic results of super- 
symmetric quantum mechanics are needed. We restrict our attention to systems with 
a quantum- mechanical hamiltonian Hs and two supersymmetry charges Q and Q\ 
obeying the defining relations of the superalgebra sqm(2): 

{Q,Q^} = Hs , [Q,Hs] = [Q\Hs] = . (16) 

The representation of this algebra relevant here is two dimensional and may be 
parametrized as |Tl|, |1^ 



There are two component hamiltonians in this system and two associated Hilbert 
spaces. If the bosonic hamiltonian acts on wave functions while the fermionic 
hamiltonian H~ acts on then the corresponding Schrodinger equations can be 
written 



^± = E±^± , (1^ 



where A is the operator A = —idy — iU' and where the supersymmetric partner 
potentials are defined by V^{y) = U'^ =F U", with U' = dyU{y) a specified function 
called the superpotential. 

The ground state of a supersymmetric system lies in the bosonic sector and has 
zero energy. Every state in the bosonic sector other than the ground state is de- 
generate with a distinct state in the fermionic sector, and the operators Q, map 
between these paired states. 

HE. Supersymmetric Coulomb System 

To construct the supersymmetric Coulomb system, the bosonic-sector combination 
— from Eq. ([18|) is identified with the radial equation (^. In the latter, a 
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suitable constant must be added to the energy eigenvalues and incorporated in the 
potential to ensure that the ground-state energy is zero. Thus, the eigenvalues are 

" 4(/ + 1 + 7)2 ~ 4(n + 7)2 ' ^^^^ 
and the bosonic potential function is 

^ ' ?/ 4(/+ 1 + 7)2 ^ ' 



The last term in Eq. ( poD is the energy shift ensuring a zero ground-state energy. 
Since it must be constant, / must be fixed to define a supersymmetric partner. 



The superpotential is specified by the function 113 



u(y) = —--i ^_(/ + ^ + i)ln2/ . (21) 

2(/ + 7 + l) ^ ' ^ ^ ^ ' 

The fermionic hamiltonian and hence the associated fermionic radial equation can 
then be calculated as 

f (f (/' + 7)(/' + 7 + 1) 1 1 1 , , 

J +^ '^^^ L-- + ———\w-{y) = E-,,w-{y) , (22) 



I dy"^ y2 y 4(//_^^)2 

where V = I + 1 and n' takes on all values of n except the lowest one, / + 1. The 
fermionic wave functions w~{y) have the same functional form as the bosonic wave 
functions W^^y) = w{y) given in Eq. (^), but with n and / replaced by n' and /' : 

Wd,n',l'iy) =Wd,n',l'{y) ■ (23) 

The two sets of eigenvalues are degenerate for n' = n: E~,^^ = E^^ . 

With fixed /, the bosonic stack of eigenstates in order of increasing energy con- 
sists of the series of kets |?t, = / + 1, I), \n = I + 2, I), . . ., with lowest energy zero. 
The associated fermionic stack has angular momentum greater by one unit, and 
starts with lowest energy corresponding to the second state of the bosonic sector: 
\n' = I + 2, / + 1), |n' = / + 3, / + 1), . . .. The use of n' here is consistent with spec- 
troscopic notation. For example, when the s orbitals of lithium are interpreted as 
the supersymmetric partner of the hydrogen atom, n' = 2 corresponds to the ground 
state, as expected. 
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A useful one-to-one correspondence between these two stacks identifies the lowest 
states with each other, and successively higher states of the bosonic sector with succes- 
sively higher states of the fermionic sector. It is defined by the following replacements 
in the bosonic wave function: 

n I — > n' = n + 1 , 

/ , — , I' = 1 + 1 . (24) 

This stack correspondence relates eigenstates with different eigenvalues. Along with 
similar stack correspondences defined below, it plays a useful role in the analyses to 
follow. 

IIF. Supersymmetric Oscillator and Composition Mapping 

The bosonic component of the supersymmetric oscillator can be obtained from 
the radial equation (|]) under a suitable energy shift. The eigenvalues are 

F+^ = 2(iV-L) . (25) 

The superpotential is specified via the function 

U{Y) = \Y^ -{L + T + l)\nY , (26) 

which generates the fermionic equation 

{-^ + ^"^^^^^^yl^^^'^ +Y^- (2L' + 2r - 1)} W-{Y) = F-„,W-iY) , 

(27) 

where L' is defined by L' = L + 1. The principal quantum number A^' takes the 
values A^' = L', L' + 2, L' + 4, .... The fermionic wave functions W~{y) have the 
same functional form as the bosonic wave functions W^{Y) = W{Y) of Eq. (P), with 
A^, L replaced by A^', L', respectively: 

WB,N',L'iY)=WD,N',L'{Y) . (28) 

The fermionic energies F^,^, = 2{N' — L' + 2) are degenerate with the bosonic energies 
for N' + l = N,F^>^,=F+^^,^,r. 



9 



In order of increasing energy, the bosonic stack with fixed L consists of the kets 
\N = L, L), \N = L + 2, L), \N = L + 4, L), . . ., with lowest energy zero. The as- 
sociated fermionic stack has angular momentum one unit greater and contains the 
kets \N' = L' = L + 1, L + 1), |A^' = L + 3, L + 1), . . .. 

The D = 1 case is unusual and warrants special attention. As mentioned above, 
the 'angular momentum' L for the one- dimensional oscillator takes the values zero 
and one, corresponding to even and odd parity. The system resembles a single stack, 
but is composed of two interlocking substacks. As a result, the spacing between 
neighboring eigenvalues is half its value in higher dimensions. Also, in constructing 
the supersymmetric partner, the energy shift is L dependent. Consequently, distinct 
shifts appear for each substack. The formalism thus establishes two independent 
supersymmetries, each of which respects the parity and only one of which may be 
considered at a time. These supersymmetries for the one-dimensional radial oscillator 
differ from the usual one for the full one-dimensional oscillator, where a single energy 
shift is effected and states of opposite parity are degenerate under the supersymmetry. 

In later sections, for reasons that emerge from the construction of the generalized 
mapping, it is more natural to focus on the supersymmetric partner of the fermionic 
oscillator rather than the fermionic oscillator itself. This system, which we call 'second 
fermionic,' has wave functions W^j^„yi given by 

WBN"L"{y) = WD,N",L"{y) ■ (29) 

Here, L" is defined by L" = L' + 1, and A^" takes values L", L" + 2,L" + 4, . . .. The 
differential equation for this system has the same functional form as the fermionic 
equation (|27|), except for the replacement of A^ and L with A^" and L", respectively. 

The oscillator bosonic sector may be put into one-to-one correspondence with the 
second-fermionic sector by making the following replacements in the bosonic wave 
function: 

N I — > N" = N + 2 , 

L I — > L" = L + 2 . (30) 
By composition of this mapping and the ones given in sections IIC and HE, a 



10 



correspondence may be established between the fermionic sector of the Coulomb 
system and the second-fermionic sector of the oscillator. It is given by 



= , / , , (32) 

(n' + 7) ' ^ ^ 

A^" = 2n'-2 + X , (33) 

L" = 21' + X . (34) 



The dimensions are still related as in Eq. (|Tj). See Figure 1. We emphasize that 
this commutative diagram involves mappings different from those presented in ref. 
IP, where the second-fermionic sector is not considered. 

III. Generalized Supersymmetry-Based Quantum-Defect Theory 

In this section, we introduce analytical SQDT for the Coulomb and oscillator sys- 
tems in arbitrary dimensions. From the present perspective, the goal is to obtain 
effective radial equations that offer sufficient flexibility to obviate the dimension and 
angular-momentum constraints of the usual mappings, while maintaining eigensolu- 
tions with analytical structure comparable in simplicity to those of the Coulomb and 
oscillator systems. 

The existence of suitable deformations of the Coulomb and oscillator systems 
satisfying these criteria is by no means apparent a priori. In what follows, we take as a 
guide the SQDT that is known to provide a useful analytical description of the valence 



structure of physical atoms in terms of an effective one-particle radial equation [|I4 
This model determines an effective radial potential modifying the three-dimensional 
radial Coulomb equation. It generates solutions with physical eigenvalues given by 
the Rydberg expression 

= -Eo/An*^ . (35) 

Here, n* is the principal quantum number modified by subtracting the quantum defect 
6, which in general depends on the angular momentum and the principal quantum 
number. In section IIIA, we generalize this model to the li-dimensional situation. For 
simplicity, we take 6 and its generalization in arbitrary dimensions to be independent 
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of the principal quantum number. This approximation is excellent in, for example, 
real alkali- metal atoms p2| . 

Similar ideas can be implemented for the radial equation of the D-dimensional 
oscillator. The resulting oscillator SQDT are presented in section IIIB. A possible 
physical application of these oscillator models is to the description of a valence particle 
in geonium atoms formed by a group of charged particles bound in a Penning trap 



IIIA. Generalized SQDT for the Coulomb System 

Given the (i-dimensional Coulomb radial equation (^) with fixed angular momen- 
tum /, we seek to implement two modifications via an effective potential Vcs{y) added 
to the left-hand side. The first desired modification is a shift in dimension, from d 
to d* = d + j, where j is an integer that in principle could depend on /. We require 
d* > 1, so j must satisfy j > 1 — d. The second desired modification is a shift in en- 
ergy eigenvalues from En^^ in Eq. (Q) to the (i*-dimensional extension of the Rydberg 
series (see Eq. (^) below). We want both these changes to be implemented while 
maintaining analytical eigenfunctions with form similar to those in Eq. @. 

Remarkably, these goals can be accomplished with a relatively simple effective 
potential, given by 

_ {n + ^)'-{n*+Y? ^ (/*+7*)(/* + r + l)-(^ + 7)(/ + 7 + l) ...^ 

Here, the quantity 7* is defined by 7* = {d* — 3)/2. The quantities n* and /* are 
defined as 

n* = Us — S = n + i — S , (37) 

/* = l + z~6 , (38) 

where 5 is the quantum defect determining the energy shifts for the generalized Ryd- 
berg series and where i = is an integral- valued function of the angular momentum. 
In the supersymmetric interpretation for the valence electron of physical atoms, is 
the number of filled lower levels with angular momentum /. The introduction of is 
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motivated by the three-dimensional case, where it is equal to the principal quantum 
number and takes conventional values in the standard spectroscopic notation. It sat- 
isfies Ug = n + i, where n takes the usual values characteristic of the exact Coulomb 
system. As an example, the s states of the supersymmetric sodium atom in three 
dimensions have i{0) = 2, giving Ug = 3, = 4, ns = 5 for the first three levels [p!3| . 
The corresponding values of n are n = 1, n = 2, n = 3. For the supersymmetric 
partner of the exact Coulomb system, i = 1 and so rig = n + 1 = n' , consistent with 
our previous notation for the supersymmetric case. 

The first term of Eq. (^) has the effect of shifting the energy levels, while the sec- 
ond term performs a corresponding shift in the angular-momentum barrier. The com- 
bined effect of both terms incorporates the desired dimensional shift. With a nonzero 
quantum defect 6{l), the effective potential fefr(z/) plays the role of a supersymmetry- 
breaking potential. The resulting radial equation has analytical solutions given in 
terms of the usual Coulomb solutions WdM,i{y) by Wd* ^n*,i*{y)- These solutions exist 
for > Z + 1, or TT-s > / + i + 1. Requiring the existence and orthonormalizability of 
the wave functions restricts 6 — i according to 

(5 - z < / + 7 + 1 + ij . (39) 

It is convenient to define a dimensionless quantity a{l) by 

a{l) = i-5+y . (40) 

The eigenvalues of the differential equation can then be expressed as 

Eq 4(n* + 7*)2 4(n + 7 + a)2 

In this equation, we have chosen the eigenenergies so that the limiting case with d = 3 
and i = j = reproduces the Rydberg series (^). For i, 6, and j chosen so that 
a = 0, we obtain the bosonic equation of the Coulomb system discussed in section 
HE, up to an energy shift. If a = 1, the fermionic sector of the Coulomb problem 
is generated instead. Moreover, the supersymmetric partner of the fermionic sector 
is generated by setting a = 2, and each successive iteration of the supersymmetry 
increments a by one unit. 

13 
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A useful stack correspondence can be established between the spectrum of the 
bosonic sector of the supersymmetric Coulomb system and the SQDT Coulomb spec- 
trum. The map is given by making the following replacements in w'^^f 

d I — > d* = d + j , 

n I — > n* = n + i — 6 = ns — 6 , 

I , — , I* = i + i-s . (42) 
IIIB. Generalized SQDT for the Oscillator System 

The techniques of section IIIA can also be applied to the radial oscillator system 
in D dimensions. For fixed angular momentum L, we can obtain an effective potential 
Veff{Y) to be added to Eq. (||) that maintains analytical eigenf unctions while inducing 
an integral shift to a new dimension D* = D -\- J > 1 and simultaneously modifying 
the oscillator energy eigenvalues via a shift to a new principal quantum number A^*. 
We refer to the resulting theory as the oscillator SQDT. 

The appropriate choice of effective potential is 

VMY) = 2(iV-iVtf-r-)+ + + + ;> - + + " + ^) , (43) 

where F* = {D* — 3)/2 and the shifted quantum numbers are given by 

= Ns-I-A = N + I-A , (44) 
L* = L + I-A . (45) 

Here, / = I{L) is an integral-valued function, analogous to in the Coulomb case, 
that can be interpreted as the number of inaccessible lower levels. The quantity 
A(A^, L) is the oscillator equivalent of the Rydberg quantum defect 6{n, I), modifying 
the radial-repulsion term in the differential equation. For simplicity in what follows, 
we take A to depend only on L, thereby paralleling the case of alkali-metal atoms for 
which S depends only on /. We have also defined a quantity playing the role of the 
principal quantum number in the spectroscopic notation, given by Ng = N+2I. If the 
dimension is unmodified and A = 0, the choice 1 = 1 yields the fermionic sector of the 
supersymmetry discussed in section IIF. In this limit A^' = A^* = N + 1 ^ Ng, which 
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differs from tlie supersymmetric limit of tlie Coulomb SQDT wliere n' = n* = Tig. 
Witli our definitions, degenerate levels in the bosonic and fermionic sectors have 
values of differing by 2/ units, but have the same value of Ng. 

The first term in Eq. (J3) implements the eigenenergy shift to the oscillator ana- 



logue of the Rydberg series, while the second term is the corresponding anharmonic 
modification to the potential that maintains analytical eigensolutions. The eigenfunc- 
tions solving the resulting effective radial equation are given in terms of the oscillator 
wave functions Wi:,N,L,iX) of Eq- ® by W^z3*,Ar*,L*(^). The existence of these so- 
lutions requires that the principal quantum number takes the values A^^ = + 2/, 
L + 2J + 2, L + 2/ + 4, . . ., or = L, L + 2, L + 4, . . .. Requiring orthonormalizability 
of the wave functions restricts the range of A — / to 

A-J<L + r + | + iJ . (46) 

We can again introduce a useful dimensionless quantity A{L) by 

A{L) = I-/\ + \J . (47) 

The eigenvalues of the differential equation can be expressed as 

F AT* P* 

' = 2Ar* + 2r* + 2A + 3 = 2A^ + 2r + 4A + 3 . (48) 

We have chosen the ground-state eigenenergy in analogy with the Coulomb case (|4l|) . 
The extra factor of 2A appears to ensure that the bosonic and fermionic spectra of the 
limiting supersymmetric case with A = 0, J = have the characteristic degenerate 
pairing. If J, A, and J are selected so that A = 0, then this SQDT system reduces 
to the bosonic oscillator discussed in section IIF, up to an energy shift. If A = 1, 
it reduces instead to the fermionic partner. If A = 2, the second-fermionic sector of 
the supersymmetric oscillator is produced. Each further iteration of supersymmetry 
produces an additional unit increment of A. Note that for fixed L the spacing between 
successive eigenvalues is always four units, regardless of the value of A. 

A correspondence can be established between the oscillator SQDT and the bosonic 
sector of the supersymmetric oscillator. The images of the wave functions are obtained 
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by making the replacements: 

D I — > D* = D + J , 

N I — . A^*=A^ + /- A , 

L I — > L* =L + I-A . (49) 

IV. Mappings between Bound States of the Coulomb and Oscillator SQDT 

Composition of the mappings in sections IIC, IIIA and IIIB allows us to establish a 
correspondence between the (i*-dimensional Coulomb SQDT and the D*-dimensional 
oscillator SQDT. This mapping is described in section IVA. One of its striking features 
is that the odd-dimensional oscillator can be imaged. In section IVB, we illustrate 
the mapping with examples involving the three-dimensional Coulomb and oscillator 
systems. 

IVA. The General Case 

The general mapping is given by 

Wn^,N.,L4Y) = K,.^^,^,_.j^^Y-'^w,.,n*,i^{{n* + ^*)Y') , (50) 

= y/in* + ^*) , (51) 

D* = 2d* -2 - 2X + J -2j , (52) 

N* = 2n*-2 + \-lJ + j , (53) 

L* = 2r + X-\j + j , (54) 

A = 2a . (55) 

This mapping, like the Coulomb-oscillator case discussed in section IIC, is based on a 
quadratic relationship between the radial variables of the two systems. The constant 
K can be chosen to preserve the normalization of the wave functions, in which case 
it has the functional form given in Eq. (plSj). Note that Eq. (p3D is equivalent to a 
generalization of Eq. (|I5D , given by 



I +4a . (56) 



-E . * 
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Note also that the allowed ranges of the quantum defects given in Eq. (|39D and Eq. 
(^) are compatible with Eq. (0), which guarantees that the image of any orthonor- 
malizable Coulomb radial system is an orthonormalizable oscillator. 

To gain insight about the flexibility of this mapping, consider the choice j = 0. 
Then, D* lies in the range l + J<D*<2d — 2 + J with allowed values separated 
by two units. Since we require J > 1 — D, any D* > 1 is possible. Moreover, d may 
take any value greater than one. The above general mapping therefore relates any 
Coulomb dimension d > 1 to any oscillator dimension D* > 1. In particular, Eq. ( ^21) 
shows that D* is odd if J is chosen to be an odd integer. This is in striking contrast 
to the usual restriction of D to even values only, as given by Eq. ([T^ . 

The Coulomb-oscillator mappings defined earlier are special cases of our general 
mapping. The correspondence of section IIC is recovered by setting i = 6 = j = 0m 
the Coulomb system and / = A = J = Oin the oscillator system, so that A = 2a = 0. 
Figure 2 is a commutative diagram showing the relationship between this simpler 
mapping and the general mapping. Similarly, the mapping of section IIP between the 
fermionic sector of the supersymmetric Coulomb system and the second-fermionic 
sector of the supersymmetric oscillator is reproduced with the choices i = 1, 6 = 0, 
j = and / = 2, A = 0, J = 0, so that A = 2a = 2. 

For the case of constant nonnegative integral A and a, Eq. ( |55| ) controls the rela- 
tionship between the supersymmetric sectors of the two systems. While any iteration 
of the supersymmetry for the Coulomb system can be taken, only even iterations of 
the oscillator supersymmetry appear. This is why we introduced the second-fermionic 
sector of the supersymmetric oscillator in section IIF. It is therefore possible to com- 
bine Figures 1 and 2 in a single commutative diagram. Moreover, the general mapping 
shows that Figure 1 can be extended downward to incorporate higher iterations of 
the supersymmetry. The result is an infinite series of mappings relating Coulomb 
systems with a = 2, 3, 4, ... to oscillator systems with A = 4, 6, 8, . . ., respectively. 

IVB. Three-Dimensional Coulomb and Oscillator Systems 

To obtain further insight about the content of the general mapping of section 
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IVA, we next restrict attention to the special case where both systems are three- 
dimensional. Since this choice can be implemented with j = 0, we assume this in 
what follows. 

The general mapping becomes 

W^,N,^L'iY) = K,^^,^^Y-'^ws,n',i*{n*Y') , (57) 

N* = 2n*-| , (58) 

L* = 21* + 1 , (59) 

A- 1 = 2{6-i) + X-l . (60) 

The values of A allowed by Eq. (^) are A = 0, 1. The orthonormality requirements 
Eq. ( ^9|) and Eq. (^) become 

S-i < l + l , (61) 

A- 1 < L + l . (62) 

We can regard Eq. (|6TD and Eq. (|6^) as conditions limiting the choice of quantum 
defects in the two systems to a semi-infinite region of the {A — I) versus {S — i) plane. 
The condition (|60D then further restricts the choice to a straight line in this region. 

One interesting special case is obtained by requiring that the oscillator be exact 
in the sense that A — / = 0. Then, Eq. (|60|) becomes 

= -a = l(i-A) , (63) 

showing that a nonzero defect in the Coulomb system is necessary. The eigenvalues 
of the equations are 

FA,,r = F„(2iV + 2A + 2) , (65) 

and the relationships among the principal quantum numbers and the angular mo- 
menta become 

L = 2l* + l , (66) 
iV = 2n*-| . (67) 
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Selecting A = 1 for definiteness, we see that Eq. (|66|) maps each successive Coulomb 
angular momentum /* = j, I* = |, /* = |, ... to every second oscillator angular 
momentum starting at L = A: L = 1, L = 3, L = 5, . . .. The mapping there- 
fore preserves the degeneracy of states. For instance, the kets \n* = |, /* = |) and 
\n* = |, ^* = |), which are degenerate in the Coulomb system, are mapped to the 
degenerate states = 3, L = 1) and = 3, L = 3) in the oscillator system. This 
feature is also a characteristic of the original mapping of section IIC. The main dif- 
ferences here are that the Coulomb effective angular momenta /* are nonintegral and, 
more importantly, that both systems are three dimensional. 

A second case of interest is obtained when the Coulomb system is exact, i.e., 
6 — i = 0. The condition j = implies that a = too. Then, Eq. (|55|) becomes 
A — I = A — ^, showing that a nonzero defect is again needed, this time in the 
oscillator system. The eigenvalues are 

^".7 = ' (68) 

F^.p. = Fo(2iV-2A + 4) , (69) 

and the mapping gives 

L* = 21 + 1 ' C^O) 

N* = 2n-| . (71) 

The first of these equations shows every second oscillator angular momentum is im- 
aged, which again preserves the degeneracy of states. If, for example, we choose 
A = 0, then the degenerate Coulomb states |r;, = 3, Z = 1) and \n = 3, I = 2) map to 
the degenerate oscillator states |A^* = |, L* = |) and |A^* = |, L* = |). 

In the above examples, the quantities A — J and 6 — i are constant. In physical 
systems such as alkali-metal atoms, 6 — i depends on I and tends towards zero as / 
increases. This feature can also be incorporated in our general mapping. It implies 
a dependence of A — / on L, which might reflect a realistic feature of a physical 
oscillator such as a cloud of particles caught in a Penning trap. 

As an example, we map the SQDT radial equation for the physical sodium atom 
into a three-dimensional SQDT oscillator. In sodium, the inaccessibility of the levels 
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occupied by the ten inner electrons is implemented by the choices i{0) = 2, = 1, 
and i{l > 2) = 0. The quantum defects 6{l) in this case are known Choosing for 
definiteness A = 1 and selecting /(O) = 2, /(I) = 1, I{L > 2) = 0, the values of A 
can be found from Eq. (|60D. Table 1 lists the results. As expected, the values of A 



tend towards ^ as L increases. 

V. Mappings for Continuum States 

In previous sections, we have explored mappings between the bound states of the 
Coulomb and oscillator systems. It is natural to consider whether similar mappings 
exist taking the unbound Coulomb states into an appropriate oscillator. This question 
is of lesser physical interest at present, so we restrict ourselves here to a brief outline 
of a possible approach to this issue. 

The Coulomb problem with energies E > can be viewed as a scattering problem. 
Following the general procedure of section IIA again yields the differential equation 
but with E > 0. In terms of the confluent hypergeometric function iFi, the 
solutions are 



Wd,E,i{y) ^ y^^'^^^exp(^±iy^j^^ 



X iFi(/ + 7 + l±^y^,2(/ + 7 + l),T2^|/^-^l . (72) 

The upper and lower signs correspond to outgoing and incoming waves, respectively. 
The results of section IIA can be recovered by taking E to be negative and choosing 
the upper sign. 



It turns out that the appropriate image oscillator system is the inverted os- 
cillator, with potential U{R) = —^MQR'^. This system is unbound. The procedure 
of section IIB gives a differential equation identical to Eq. (^ except that the sign of 
the potential is reversed. The solutions W£)^f,l(X) are 



WD,F,L{y) ^ exp {±ltY^) iFi (|(L + r + I) T L + r + |, TtY' 

(73) 

These functions may not be physically permissible, but are relevant for the purposes 
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of establishing a mapping. With the choice of the upper sign, the wave functions for 
the usual oscillator may be obtained up to a constant by the analytic continuation 
iY^ and F -iF. 

A correspondence analogous to the mapping of section IIC exists between the 
continuum Coulomb states and the inverted-oscillator functions. It is 

Wn,F,L{Y) « Wd,E,i I I , (74) 



y2 = 2y^E/Eo , (75) 
D = 2d-2-2X , (76) 
L = 21 + X , (77) 

F / En 

There are many similarities between this mapping and the one discussed in section 
lie. Again, A must be integral, so only even- dimensional oscillators are available as 
images. Also, the angular momenta L are restricted to being either all odd or all even, 
thus eliminating half the oscillator states. However, the energy relation (^), unlike 
(p!5|), involves continuous values of E and F. It also reveals that negative energies 
F are excluded from the mapping. As one energy tends to zero the other tends to 
infinity. 

As an aside, we remark that the negative energies F do appear when considering 
the repulsive Coulomb problem. The differential equation of this problem is mapped 
into the inverted-oscillator differential equation |24| by a map with (^) through (77) 
unchanged but with a negative sign taken for the square root in Eq. (ff^). 

Although this lies outside the scope of the present work, it seems feasible that 
supersymmetry could be introduced into these systems along with the corresponding 
SQDT. We conjecture that this allows for odd dimensions D. Since parabolic coor- 
dinates have some advantages for scattering problems, it would also be interesting to 
perform an analysis in terms of the dual parabolic-coordinate supersymmetries of ref. 
instead of the spherical-coordinate supersymmetry used here. 
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VI. Summary 

In this paper, we generalized the radial mappings first identified by Schrodinger 
that relate the Coulomb and oscillator systems. Our principal result is a mapping 
between the supersymmetry-based quantum-defect theories for the Coulomb and os- 
cillator systems in arbitrary dimensions. In particular, odd oscillator dimensions can 
be accessed as well as the usual even ones. The mapping and some of its limits are 
illustrated in Figures 1 and 2. 

In deriving this result, we have extended to arbitrary dimensions the analytical 
SQDT in three dimensions used to describe physical alkali-metal atoms. An analogous 
SQDT for the harmonic oscillator in arbitrary dimensions has also been presented. 
In suitable limits, these theories reproduce the bosonic and fermionic sectors of the 
corresponding supersymmetric quantum-mechanical systems. We have elucidated a 
basic relationship between the supersymmetric radial Coulomb and oscillator systems: 
the gth iteration of supersymmetry for the Coulomb system corresponds naturally to 
the 2qth iteration of supersymmetry for the oscillator. For the special case of the one- 
dimensional radial oscillator, we uncovered a quantum-mechanical supersymmetry 
in which the parity is restricted to be cither odd or even. We have also briefly 
considered mappings relating the continuum-spectrum states of the Coulomb and 
oscillator systems. 

The issue of the physical relevance of our results has also been addressed in part. 
The three-dimensional Coulomb SQDT is known to provide a good analytical descrip- 
tion of the behaviour of Rydberg atoms. Our mapping provides a means of obtaining 
an equivalent analytical oscillator SQDT. An explicit example mapping the sodium 
atom to an oscillator SQDT is given in Table 1. It is also possible that an oscillator 
SQDT could be used to describe a suitable physical system, perhaps the Penning 
trap. If this can be realized in practice, the generalized mapping presented here could 
provide a connection between two apparently disparate physical systems. 
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Table 1. Possible parameters for a mapping between the radial Coulomb and oscillator 
systems in three dimensions. For the choices j = and A = 1, values of /, i, n, Ug, 
and 6 are tabulated for sodium along with the corresponding values of L, I, N, Ng, 



and A under the mapping (|60D. The quantities I{L) have been selected to fill all 



levels below Ng = 5. 



Figure 1. Supersymmetric mappings. Relationships are shown interconnecting the 
bosonic and fermionic partners of the Coulomb and oscillator systems. The diagram 
is commutative. 



Figure 2. SQDT mappings. Relationships are shown interconnecting the bosonic and 
SQDT sectors of the Coulomb and oscillator systems. The diagram is commutative. 
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Coulomb system (sodium) 


Oscillator system 




i 


n 


Us 


5 


L 


/ 


N 




A 





2 


> 1 


> 3 


1.35 


1 


2 


> 1 


> 5 


1.20 


1 


1 


> 2 


> 3 


0.859 


3 


1 


> 3 


> 5 


1.218 


2 





> 3 


> 3 


0.01 


5 





> 5 


> 5 


0.52 


3 





> 4 


> 4 


0.00 


7 





> 7 


> 7 


0.50 


> 4 





> / + 1 


>/ + l 





> 9 





> L 


> L 


0.5 



Table 1 
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Coulomb 
bosonic 

\n,l) 



2/ + A = L 



2n - 2 + A = TV 



Oscillator 
bosonic 
\N,L) 



I' = 1 + 1 
n' = n + 1 



L" = L + 2 
N" = N + 2 



Coulomb 
fermionic 

W, I') 



2/' + A = L' 



2n' - 2 + A = A^' 



Oscillator 
2nd fermionic 
\N",L") 



Figure 1 
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Coulomb 2/ + A = L Oscillator 

bosonic — ^ bosonic 

I 7\ 2n- 2 + A = iV ^\ 

|n,0 2d-2-2X = D 1^'^) 



/* = / + z- J L* = L + /- A 

n* =ns-S N* = Ns-I-A 



2l* + X + j -y = L 



Coulomb 2n* - 2 + A + 7 - i J = N* Oscillator 

2d'-2-2A-2, + J-g' ^'^^ 
\n',l') \N\L') 



Figure 2 



28 



